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ABSTRACT 

Context. Current solar energetic particle (SEP) propagation models describe the effects of interplanetary plasma turbulence on SEPs 
as diffusion, using a Fokker-Planck (EP) equation. However, EP models cannot explain the observed fast access of SEPs across the 
average magnetic field to regions that are widely separated in longitude within the heliosphere without using unrealistically strong 
cross-field dilTusion. 

Aims. We study whether the recently suggested early non-diffusive phase of SEP propagation can explain the wide SEP events with 
realistic particle transport parameters. 

Methods. We used a novel model that accounts for the SEP propagation along field lines that meander as a result of plasma turbulence. 
Such a non-dilfusive propagation mode has been shown to dominate the SEP cross-field propagation early in the SEP event history. 
We compare the new model to the traditional approach, and to SEP observations. 

Results. Using the new model, we reproduce the observed longitudinal extent of SEP peak fluxes that are characterised by a Gaussian 
profile with <r = 30 - 50°, while current diffusion theory can only explain extents of 11° with realistic diffusion coefficients. Our 
model also reproduces the timing of SEP arrival at distant longitudes, which cannot be explained using the diffusion model. 
Conclusions. The early onsef of SEPs over a wide range of longitudes can be understood as a result of the effects of magnetic field¬ 
line random walk in the interplanetary medium and requires an SEP transport model that properly describes the non-diffusive early 
phase of SEP cross-field propagation. 
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1. Introduction 

Solar energetic particles (SEPs) are accelerated up to relativistic 
energies during solar eruptions. Their propagation in interplane¬ 
tary space is controlled by the large-scale Parker spiral magnetic 
field. SEPs can also propagate across the mean field as a result of 
large-scale drifts ( Marsh et al.| 2013| l and the turbulent magnetic 
field fluctuations superposed on the mean field. 

SEP propagation in interplanetary space is typically mod¬ 
elled by solving a Eokker-Planck equation that describes the 
effect of the plasma turbulence on the particles as diffusion 
( [Parker |1965| jJokipiij |1966| l. The field-aligned propagation is 
affected by small-scale inhomogeneities ( [Jokipii |1966| ), whereas 
the cross-field propagation is described as diffusion caused by 
random-walking field lines Pokipiij |1966| [Matthaeus et al.j 
2003 Shalchi|[2010| l. These theories give the ratio of the cross¬ 
field and field-aligned diffusion coefficients ~ 0.01 near 

Earth, and these values are supported by cosmic ray observa¬ 
tions ( [Burger et al. [2000 [ Potgieteretal.j 20141 and particle sim- 
ulations ( [Giacalone & Joki^ 1999| l. Recently, several works 
have studied the effect of cross-field diffusion on SEP evolution 


([Zhang et al. 2009[[Dr6ge et al. 2010[ He et al.[ 201 1[ Giacalone 
& Jokipiij 2012 Qin et al.J 2013| l and showed up to 360° helio- 
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longitudinal extents for wide source regions with onsets within 
days of the SEP injection. 

Recent multi-spacecraft SEP observations, made with the 
near-Earth spacecraft SOHO and ACE and the two STEREO 
spacecraft, show that SEPs from a single solar eruption have fast 
access to a very wide range of longitudes, even 180° from the 
eruption location (e.g. Dresing et al.[ 2012[ Richardson et al.j 


20141. Similar event extents as measured by longitudinal peak 


intensity distribution width have been observed both in gradual 
and impulsive events (Wiedenbeck et al. |2013[ Cohen et al.j 
[2014[ [Richardson et al.[ [2014] l, which in traditional view have 
very different acceleration origins (e.g. Reames 1999| l. Detailed 
analysis of two SEP multi-spacecraft-observed events, with SEP 
transport fitted with a 3D EP transport model, favour a small 
source region of about 20° for the SEPs and a high diffusion co- 


efficient ratio a'j./a'h ~ 0.1 (Dresing et al. 

2012 

Droge et al. 

20141. High ratios in the range of 0.1-1 have also 

been inferred 

in other studies of heliospheric particles (e.g. [Zhang et al.| 2003 


Dwyer et al.[ 1991) . 


High values of Kj_Ik\\ are not supported by our current the¬ 
oretical understanding of SEP transport in plasma turbulence, 
however. At 1 AU, the amplitude of transverse magnetic fluctu¬ 
ations, normalised with the background magnetic field, 
is observed to be between 0.1 and 0.5 (e.g. Burlaga & Turner j 


19761, for which full-orbit particle simulations suggest Kj_Ik\\ < 
0.01 [Giacalone & Jokipii|([T999|l. Current theoretical understand- 


1 











































































































T. Laitinen et al.: Solar energetic particle access to distant longitudes through turbulent field-line meandering 


ing also supports much lower values for the ratio (e.g.|Pei et al.l 

[2^ . 

In this paper, we introduce a new model for heliospheric 
SEP propagation that is capable of reproducing the fast trans¬ 
port in longitude that is observed at the beginning of SEP events. 
Our work is based on the notion that SEPs spread initially much 
faster than allowed by a diffusion description, and that this is a 
result of field-line meandering ( [Laitinen et al. 2013| l. The cross- 
field propagation of SEPs along meandering fields has been 
analysed previously to study SEP intensity dropouts (Giacalone 


et al. 


20001 and to determine diffusion coefficients (Kelly et al. 


2012 I. Laitinen et al. (|2013jl showed for a uniform magnetic field 


that particle propagation across the mean field is initially a fast 
non-diffusive process along the meandering field-lines that only 
later relaxes to diffusion. Our new approach combines field-line 
meandering and time-asymptotic diffusion in a new description 
of early-time SEP propagation in the Parker field, the EP-i-PLRW 
model. We compare the new model to the traditional diffusion 
approach (the EP model). 

2. Models 

2.7. FP model 

In the traditional approach, SEP propagation is solved using the 


Eokker-Planck equation for the distribution function / (Roelof 


[T969i [Skilling [TWTl jlsenberg] [T997i [Zhang et al.j [2009 i 
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where v, p and fi are the particle speed, momentum, and pitch- 
angle cosine, and b the solar wind velocity and a unit 
vector along the local mean magnetic field, respectively, and 
L - -BKdBjds) the focusing length of the particles, with s the 
arc length along the field line, and Q is the injection function. 
The terms of the order (dVswldt)fv ( [Isenberg \991\ are typi¬ 
cally neglected in SEP transport studies as small compared to 
the other terms (e.g. Zhang et al. 2009 [ Droge et al. |2010| l, and 
they vanish for the constant radial solar wind considered in our 
study (see Sect.[2.3[). 


The pitch-angle diffusion coefficient and the cross-field 


diffusion coefficient are given as ([Jokipii[ [1966[ [Matthaeus 

[er^|2^ 
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The field-aligned spatial diffusion tensor is given as 

' 0 0) 

0 a:j_ 0 
,0 0 0 , 


(5) 


that is, without parallel diffusion, because this is accounted for 
with D, 


A'A'- 


The changes in SEP energy during propagation in the inter¬ 
planetary space are given by the terms in the third line of Eq. ([^l. 
They include momentum changes that are due to the diverging 
solar wind, and in an approximate way, the effect of drifts on the 
particle momentum ( [Dalla et 'ak] [2015| l. However, the equation 
does not contain the drifts in its convective term (second term in 
the first line), which means that the transport picture presented 
by Eq. (0 is incomplete. The effect of latitudinal drifts on the 
particle energy was recently discussed by Dalla et al. ( 2015[ l, 
who found that in 100 hours a 10 MeV proton can lose almost 
half of its energy through drifting along the Ysw x B electric field 
and adiabatic deceleration. An improved approach, the drift ki¬ 
netic theory, introduces drifts within the convective term (e.g.|I^ 
[Roux & WebbI [200^ [Webb et al.|[2009| ). 

The energy changes that are due to drifts and the diverging 
solar wind may be significant when performing a detailed fit¬ 
ting of the SEP event onset and the decay phase (e.g. Ruffolo[ 
[1995[[Dalla et al.[[2015[ l. However, we here consider neither, but 
only the first ten hours of the event. Eor simplicity, the energy 
change term is therefore neglected in this study. Because the spa¬ 
tial transport associated with drift is small for a 10 MeV proton 
within the early phase of an SEP event (Dalla et al. 2015| l, drifts 
are neglected within the convective term. 

Particle diffusion parallel and across the mean field direc¬ 
tion, as understood by current theories, are not independent of 
each other. The cross-field diffusion coefficient in Eq. ([^ de¬ 
pends on the parallel diffusion coefficient, resulting in compound 
diffusion (e.g. [Kota & Jokipii 2000[ Qin et al.j 2()02[ Matthaeus 
[et al.[ [2003[ ). The pitch angle dependence of the cross-field dif¬ 
fusion coefficient has also been discussed recently, with sugges¬ 
tions for using oc \ p \ (e.g. Qin & Shalchi 2014[ l or oc ri , the 
particle Larmor radius ( Droge et al.[[2010] l. While the suggested 
different forms of pitch angle dependence of a:j_ influence parti¬ 
cle cross-field propagation (e.g. Strauss & Eichtner][2015[ l, it has 
not been clearly established which type of dependence best rep¬ 
resents particle behaviour in turbulent fields. We therefore use 
the conventional pitch-angle independent form of here. 

We solved Eq. jEl using stochastic differential equa¬ 
tions (SDEs, e.g. Gardiner! 2009} , which solve the equations 
using pseudo-particles. The method is briefly described in 
Appendix [A} further details are given in Kopp et ^P012[ l. We 
used enhanced pitch angle scattering across ju = 0 by using a 
pitch angle diffusion coefficient of the form 

D,^^Ko{\pr^ +h){1-p^), (6) 

suggested by Beeck & Wibberen^ ( [1986[ l, where a-q contains the 
//-independent terms in Eq. (0, q is the spectral index of the slab 
turbulence, and // = 0.1 is a parameter that enhances scattering 
between pitch angle hemispheres. The EP equation is solved in 
a Parker spiral magnetic field with the magnitude B given as 


where S (k) is the turbulence power spectrum, S y the spectrum of 
slab turbulence (see Sect. [2.3[ l, Q and are the particle Larmor 
freque ncy and radius, respect ively, the free parameter gnlcc - 
1 / V3 (Matthaeus et al. 2003|, and the parallel mean free path 


B(r) = Bo - 



(7) 


where Bq = 5 nT is the magnetic field at heliocentric distance 
ro = 1 AU, and a - VswUFlQ&mO), where - 400 km/s, 
Go = 2.8631 ■ 10“® rad/s is the solar rotation rate and 6 the 
colatitude. 
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2.2. FP+FLRW model 


Laitinen et al. ( 2013| l demonstrated that the early-time particle 
transport across the mean magnetic field is dominated by the 
particles propagating along meandering held lines. With respect 
to the mean magnetic held direction, the early-time cross-held 
propagation of the particles is therefore only weakly stochastic. 
At later times, [Laitinen et ahlpOlSl l showed that the cross-held 
extent of the particle population begins to widen, resulting in an 
asymptotic diffusive behaviour. 

In this study, we used an FPh-FLRW model that explicitly 
introduces the random walk of held lines across the Parker spi¬ 
ral magnetic held and particle propagation along these random¬ 
walking held lines. A similar model was successfully used by 
Laitinen et al. ( 2013| l for Cartesian geometry. In addition to prop¬ 
agating particles along meandering held lines, the particles are 
diffused across them with the time-asymptotic cross-held parti¬ 
cle diffusion coefhcient to facilitate the late-time widening of 
the particle population. Laitinen et al. ( 2013| l showed that such 
a model reproduces both the initial wide extent of the SEP event 
(where the FLRW dominates) and the asymptotic diffusive be¬ 
haviour at later times (where the particle diffusion dominates) 
well and that is signihcantly better than the FLRW or FP on their 
own. 

In our model, the held line random walk is facilitated by 
the turbulent interplanetary magnetic held, which is described 
by a power spectrum. A mix of 2D and slab modes in turbu¬ 
lent magnetic held gives rise to diffusive spreading of held lines 
( Matthaeus et al.|[l995j l. This allows us to describe the path of a 
particle following a meandering held line using SDEs, with 


drjr\i) = ^J2DFL(r\\)dr\lWJ_ 


(8) 


where is the displacement across the Parker held direction 
for advance dr\\ along the held, Dfl is the held-line diffusion co¬ 
efhcient, and Wj_ a Gaussian random number with zero mean and 
unit variance. The random-walking path is calculated for each 
simulated particle before the pseudo-particle trajectory is inte¬ 
grated. The pseudo-particles are propagated along and diffused 
from the meandering held lines according to Eq. ([^l, with cross- 
held diffusion given by Eq. ([^, pitch angle diffusion as given in 
Eqs. and (j^, and focusing, where we used the Parker held 
geometry when calculating the focusing length. 

The held-line diffusion coefficient was obtained as presented 
in [Matthaeus et al.j (19951, for composite slab and 2D turbu¬ 
lence. The contribution from the slab component is proportional 
to 5||(k|| = 0), which in our turbulence model (see below) van¬ 
ishes. We therefore only considered the 2D contribution to the 
held-line diffusion, giving ([Matthaeus et al. 1995|l 


Dfi. — 


7 d\^S,A^J/kl 

R2 


I 1/2 


Consistent with the SDE approach in the traditional model, 
the effect of the cross-held step length is not considered in the 
particle propagation time. This affects the accurate timing of par¬ 
ticle propagation (e.g. [Strauss & Fichtner|[20T4| l. 

2.3. Turbulence model 


To compare the EP and EPh-ELRW models, it is essential that 
D/ifj’ '‘'j- Dfl, as dehned in Eqs. 0, ([^ and (|^, are consis¬ 
tent with each other. We ensured this by denving them for a sim¬ 
ple interplanetary turbulence model. We considered only wave 
modes parallel (slab) and perpendicular (2D) to the mean mag¬ 
netic held (Gray et al. 1996| l and dehned the power spectrum 
as 

3(k) = 3(k,ro) = SMm\) + S\\(h)6(ki), (10) 

where ky and - |kx| are magnitudes of the wave number. The 
spectra are given as broken power laws 



’C||L||5By 

Lq < ky < L|| 


S\\(k\\) - ■ 

CyLybBjl 

ky > L"' 

(11) 




and 

CiLi6B\ 

< ki < 


Si(ki) = ■ 

CiLi5B\ 

ki > Ll\ 

(12) 





where Ly and are the spectral breakpoint scales with a value 
of 0.007 AU used for both scales, consistent with solar wind 
observations (e.g. & Marsch] |1995| and references therein). 

The largest scale in the turbulence model, Lo(r), is taken to be r, 
a natural choice for scaling in spherically expanding solar wind. 
The turbulence amplitude 6B^ is related to the normalisation fac¬ 
tors Cy and Cl to give ~ ^ f ‘^||,±(k)<7k. The turbulence 

amplitudes are normalised to give the parallel scattering mean 


trie distance (e.g. 

Palmer 

1982 

1 , and 6Bf^ : 6B\ = 20% : 80% 

(Bieber et al. 

1996 1 , giving SB^IB^ - 0.04 at 1 AU, in line with 


interplanetary turbulence observations (e.g. Burlaga & Turner 
fT976i) . 

The plasma turbulence evolves as the solar wind propagates 
from the Sun to the interplanetary medium (e.g. Tu & Marsch[ 


19951. The amplitude of the huctuations changes as a function 


(9) 


It should be noted that particle propagation as a beam along held 
lines that random-walk according to the held-line diffusion co¬ 
efficient in Eq. ([^ causes transport across the mean held much 
faster than diffusion with the coefhcient given by Eq. ([^. The 
latter aims to describe the cross-held spreading of particles at 
late times, when the pitch-angle distribution of the particles is 
isotropic. The two descriptions converge for a particle beam in 
turbulence with vanishing slab component, anlcc = 1 and a 
scatter-free velocity correlation {Vj(0)v2(f')) = in Eq. (3) of 
Matthaeus et al. ( 2003| l. The two descriptions of particle cross- 
held propagation are therefore consistent with each other in the 
appropriate limit. 


of radial distance as the properties of the background plasma 
change. The turbulence also evolves non-linearly, with energy 
cascading towards smaller scales. The radial and spectral evolu¬ 
tion of turbulence in the heliosphere have been studied in sev¬ 
eral works over several decades. Many of the works, however. 


considered isotropic or slab turbulence (e.g. Tu et al. 1984 Tu 


|1987[ [Zhou & Matthaeu^ |1990[ jVainio et al.| |2003] l instead of 

the more realistic anisotropic geometry, or the composite mode 
used in this and many other recent SEP transport studies. Some 
formulations exist for spectral transport separately to small par¬ 
allel and perpendicular scales (e.g. Cranmer & van Ballegooijen[ 
[2005 [ [Laitinenl |2005[ ). These works, however, have not been 
considered with respect to particle transport coefficients. The re¬ 
cent detailed heliospheric turbulence transport model by Zank 
et al. (2012|l provides the evolution of the total turbulence en 


ergy and correlation lengths instead of the turbulence spectrum. 
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Fig. 1. Left; the parallel and perpendicular scattering mean free paths for 10 MeV protons obtained with Eqs. (|^-Q and Aj_ - 3 kj_Iv. 
Right; Diffusion coefficient ratios. For the red curve, the cross-field diffusion coefficient is obtained as a:j_ = vDfi, representing a 
cross-field diffusion coefficient that is due to field-line meandering of unscattered particles. 


and the application of their model to SEP transport modelling is 
a challenging task that has so far not been attempted. 

Developing an SEP transport model that considers the spec¬ 
tral evolution of composite turbulence in the heliosphere is there¬ 
fore an undertaking that has not been considered before, and it 
warrants a separate study. In this work, we take a simpler ap¬ 
proach and only consider the radial evolution of the turbulence. 
We used the WKB description ( Richter & Olber^ 1974[ Tu et al. 


1984^ for the slab and 2D components. For simplicity, we did not 


consider the evolution of the spectral shape nor the wave refrac¬ 
tion, and we kept the wave components in their pure slab and 
2D geometries. Furthermore, we did not consider the change in 
k that is due to the varying wave phase speed because assuming 
a constant inertial-frame wave phase speed agrees well with ob¬ 
servations (Vainio et al. 2003 | l. With these simplifications, the 
turbulence evolves as 


5||,±(k||,i,r) = 5||,i(k||,i,ro) 


VrO VAr 


VrVArO 

= Si_i_(k\u^roW(r), 


^eQ 


1/2 


(13) 


where VAr is the radial component of Alfven velocity, Vr - 
Vsw.r + VAr and He the solar wind electron number density, with 
the subscript 0 denoting the values at reference distance tq. We 
used a constant solar wind velocity Vsw,r 0 and electron density 
ne(r) = rieo r^lr^. Using these, we find 


W(r) = 


Ujw.rO + VaO 
Vsw,r0 + 7^40 


(14) 


which is consistent with the observed approximately oc trend 
in interplanetary space (e.g. Bavassano et al. 1982| l. We here 
used Vsw,tO - 400 km/s and Va,r 0 - 30 km/s to represent the 
values at ro = 1 AU. 

It should be noted that our model for the field-line random 
walk differs from that of |Giacalone| ( |200T] l, where the field-line 
random walk is caused by the motion of the magnetic field foot- 


points at the Sun that are due to supergranulation. The Giacalone 
( |2001| l model provides an explanation for the source of the field¬ 
line random walk at the Sun, but it does not allow for turbulence 
evolution in interplanetary space, which further modifies the tur¬ 
bulence as the magnetic fields are convected with the solar wind 


(e.g. Bruno & Carbone 2005 and references therein). In our 
model, the field-line meandering is derived from the turbulence 
spectrum and the corresponding field line diffusion coefficient, 
with no limiting assumptions on the interplanetary evolution of 
the turbulence. 


2.4. Transport parameters 

After specifying a turbulence model, we now calculate the par¬ 
ticle transport parameters that are fully consistent with the prop¬ 
erties of the turbulence for both the FP and FP-hFLRW models. 
The transport parameters are shown in Fig. [T] for 10 MeV pro¬ 
tons. In the left panel we show the particle transport parameters 
as parallel and cross-field mean free paths, which are the same 
for both the FP and the FPh-FLRW model, and the right panel 
shows the ratio of the cross-field and parallel diffusion coeffi¬ 
cients, The values of the particle diffusion coefficients are 


broadly similar to those used in other studies (e.g. Zhang et al. 
[200^ [PeTetaLl [Mol l. 

At small heliocentric distances, our model gives high values 
for the parallel mean free path (left panel of Fig. [^l, allowing the 
particles to propagate along the field lines essentially without 
scattering from the Sun to 0.4 AU. A similar transport regime is 
also present in models that use a constant radial diffusion coeffi¬ 
cient to parametrise the parallel diffusion coefficie nt (e.g. [Zhang 
et al.[ 20091 [Drdge et al. 2010[ He et al. 2011 Giacalone & 


Jokipi [[2012[|Qin et al.[[2(H3r 


In addition to Kj_/k\\, we show also the ratio for particles that 
propagate as a beam along the meandering field line, v ■ Dpi, 
in the right panel of Fig. [T] As discussed in Sect. 2.2 v ■ Dpi 


is always larger than Kj_. In addition, as shown by (Fraschetti 
& Jokipii 201 l|l and Fraschetti & Giacalone (2012|l, for typical 


interplanetary turbulence parameters the rate of decoupling of 
particles from their field lines is always slower than the spread¬ 
ing of particles across the mean field. Thus, at early times the 
SEP propagation along random-walking field lines dominates 
over both the asymptotic perpendicular diffusion and the particle 
decoupling from the field. 

The cross-field mean free path for both particles and field 
lines is very small close to the Sun, of the order of 10“^ AU. The 
value is deceptive, however; a low value close to the Sun can 
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Fig. 2. Sample path (red curve) of a 10 MeV proton injected at the Sun (located at origin) with the FP (left panel) and FPh-FLRW 
model (right panel), respectively. The black curve shows the Parker spiral field, while the dashed blue curve shows a meandering 
field line in the FP-hFLRW model (right panel). 


result in large longitudinal spread. In radial geometry, the field¬ 
line diffusion coefficient can be written in terms of a longitudi¬ 
nal diffusion coefficient, Dfi - /(2r) = = 

r^Dfijp. In a radial geometry, a constant longitudinal diffusion 
coefficient is therefore represented by a diffusion coefficient 
that is proportional to r^. For our turbulence model, we have 
Dfl 4 > = (40°)^/AU close to the Sun, and Dfui> - (10°)^/AU at 
1 AU. 


It should be noted that a simple parametrisation with k^_I^ a 
constant or oc r does not reproduce the behaviour of the diffusion 
coefficients obtained with a full model of the turbulence, such 
as the one used in this study and in |Pei et ak ( 2010) 1. The radial 
evolution of turbulence and its effect on the diffusion coefficients 
will be discussed in more detail in a separate study. 


3. Results and discussion 

We here study a simple injection profile, 

Q{r, 6, (p, t) = 5{r - lro)(5(6» - njl) ■ 

6(n - mE - Eo)6(<pmt - to), 

for Eq - 10 MeV protons and trace the propagation of pseudo¬ 
particles through interplanetary space according to the model 
equations. The resulting profile can also be convolved with 
different injection profiles to produce response for time- and 
longitude-extended injections, which will be a subject of a fu¬ 
ture full parameter study. 

The propagation of the particles in FP and FP-hFLRW mod¬ 
els is shown in Fig. as a 2D projection of the 3D orbits. The 
red curves depict the path of a simulated pseudo-particle, and 
the black curves show the Parker spiral. In the right panel, the 
blue dashed curve depicts the diffusively meandering field line. 
In both models, the particles scatter along the field line and dif¬ 
fuse away from it in the cross-field direction. In the FP-(-FLRW 
model, however, the field line itself meanders away from the 
Parker Spiral field. As discussed above, the diffusion of the field 
line is faster than the diffusion of a particle across the field line, 
thus the particle orbit tends to follow the meandering field line 


(blue dashed curve in Fig. and as a result spreads faster across 
the mean Parker spiral field than in the FP case. 

In Fig. 1^ we compare the particle distributions predicted by 
the diffusion (FP) (left) and our new FP-hFLRW (right) model, 
one (top row) and three (bottom row) hours after injection. The 
distribution contours are calculated for particles near the ecliptic 
plane, between latitudes +10°. In the FP model, particles diffu¬ 
sively spread in longitude as a function of time, and as they are 
focused rapidly away from the Sun in the diverging magnetic 
field, the most significant spreading of the particles takes place 
at distances > 0.5 AU. In the FP+FLRW model, however, the 
particles have a much wider distribution in longitude and show 
significant spreading close to the Sun. 

To compare the models with observed SEP events, we show 
in Fig. I^the time evolution of the simulated intensities at three 
different longitudes for the two models, compared to an event 


that took place during November 3-4, 2011. Gomez-Hetrero 
et al. ( 2015| l have analysed this event using a ID transport model. 


which, based on anisotropy observations, suggests a source re¬ 
gion with 270° longitudinal extent. However, they found no sup¬ 
porting evidence for such a wide source and suggested a dis¬ 
turbed interplanetary magnetic field as one of the possible causes 
for the event extent. 

The SEP event is related to a solar eruption at longitude 50° 
east of the STEREO B (STB) spacecraft, with an M4.7 flare 
maximum at =22:41 UT on November 3 at N20E62 from 


Earth’s point of view ( Richardson et ak) 2014 Gomez-Hetrero 
|et al.| |2015| l. Using the longitude angle (p as defined in Eig. 
where the particle source is at = 0°, STB was at longitude 
<^STB = 50°, (see lower right panel of Eig. |^. In this coordi¬ 
nate system, the SEP source is connected along the Parker field 
line (solid black spiral curve in the lower right panel of Pig. 
to longitude (p - -62° at radial distance of 1 AU (red circle). 
During the event, the STEREO A (STA) and SOHO spacecraft 
were located at longitudes (psjA - -103° and 0 soho = 152°, 
respectively. 

Observed intensities of 6-10 MeV protons are shown in 
Pig.j^with dashed green curves. They were obtained with ERNE 
onboard SOHO ( Torsti et al.[ |1995| l (top left panel), LET on¬ 
board STB ( |Mewaldt et al.| |2008| ) (top right), and the LET in- 
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Fig. 3. SEP particle density in arbitrary units for particles between latitudes -10° and 10°, one hour (top row) and three hours 
(bottom row) after an impulsive injection at 1 solar radius and longitude 0 = 0°. Left panels show the result of the FP model, while 
the right panels show the FP-i-FLRW model. The black curve shows the Parker field-line starting at 0 = 0°. The field line reaches a 
radial distance of 1 AU (red circle) at 0 = -62°. 


strument onboard STA (bottom left). The blue and red curves in 
Fig. I^show the intensities from the FP and FP-i-FLRW models, 
respectively, for the injection function Q with fo = fmax - 500s, 
(subtracting 500 seconds of light-travelling time from Sun to 
1 AU). The injection intensity is fitted so that the highest in¬ 
tensities of the models match the observed SEP highest intensity 
at STA. 

The SEP event was best observed with LET/STA, located 
41° clockwise from the well-connected longitude of 0 = -62°. 
The 6-10 MeV proton intensity increased by 5 orders of magni¬ 
tude within a few hours (bottom left panel of Fig. |^. The event 
was also observed by LET/STB and ERNE/SOHO at 112° and 
146°, respectively, from the well-connected field-line. Even with 
such poor connection, SEP intensities began to rise within an 
hour of the observed SEP onset at STA, to rise by 3 orders of 
magnitude above the pre-event background. It is therefore clear 
that the SEPs had rapid access to a wide range of longitudes dur¬ 
ing this solar event. 

While we did not perform a full transport fitting in this re¬ 
port, our study implies that the rapid access of SEPs during the 
November 3^, 2011 event to the whole inner heliosphere can be 
achieved with our FP-hFLRW model even with a narrow source 
region. The FP-i-FLRW-modelled particles are observed within 
an hour at the longitudes corresponding to the November 3, 2011 
spacecraft locations. The FP model can reproduce the rapid on¬ 


set of the SEP event at the STA spacecraft, with a delay of half 
an hour compared to the FPh-FLRW model. It is not capable of 
reproducing the SOHO or STB observations, however, because 
it fails to spread particles to these longitudes to intensity levels 
above the pre-event background during the first ten hours of the 
event. 

The inability of the FP model to spread particles to the longi¬ 
tudes of the SOHO and STB spacecraft within the first ten hours 
is shown in Fig.|^ where we show the evolution of the simulated 
SEP event as it would have been measured at 1 AU at different 
heliographic longitudes, with fo = 0. For the FP model the in¬ 
tensity falls very sharply as a function of longitude, explaining 
the failure of the FP model in Fig. to describe the November 
2011 event without resorting to extremely wide source regions 

. As shown in Fig. the SEPs 
spread quickly for the FPh-FLRW model, and the SEP onset 
within a wide range of longitudes, between 0 = -150° and 50°, 
is delayed only by an hour compared to the onset at the well- 
connected longitude, 0 = -62°. Such a rapid spreading of first 
particles is a key feature of observations: the earliest SEP onsets 
at 100° from the flare longitude are delayed only an hour com¬ 
pared to onsets at best-connected longitude ([Richardson et al. 
|m4| ). 

The longitudinal extent of an SEP event is frequently de¬ 
scribed by means of a Gaussian fit, /(0) = /oexp{-(0 - 


(Gomez-Herrero et al.[ |2015 1 
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Fig. 4. SEP proton intensities as a function of time at 1 AU. The green dashed curve shows 6-lOMeV SEP proton intensities observed 
during a solar eruption on Nov 3^, 2011, at SOHO (top left panel), STEREO B (top right) and STEREO A (bottom left) spacecraft. 
The blue and red curves show the 10 MeV proton intensities simulated using the EP and EP-i-PLRW models, respectively, using 
an impulsive injection at flare maximum, to =22;41-500s UT on Nov 3 (dashed vertical line) at the flare longitude. The spacecraft 
locations are shown in the bottom right panel by filled green circle, yellow triangle, and blue square for SOHO, STEREO B, and 
STEREO A, respectively. The blue arrow denotes the SEP injection longitude, the thick solid line and red circle are the same as in 
Pig. 0 and the dashed curves show the Parker spiral field-lines connecting the spacecraft to the solar surface. 


0o)^/(2cr^)), to the peak SEP intensities observed by spacecraft 
at different longitudes. Using STEREO and near-Earth space¬ 
craft, several studies have reported values of cr = 30° - 50° to 


best represent the range of SEP event widths ( 

Lario et al. 2006 

2013 

Wiedenbeck et al.[[2013[ Dresing et al. 

2014 Richardson 

et al. 

2014[l. To compare our models to observations, we per- 


formed a Gaussian fit to the peak intensities of the simulated SEP 
events at 1 AU for the two models presented in this work. The 
longitudinal distributions and the corresponding hts are shown in 
Pig. I^with solid and dashed curves, respectively. The grey area 
depicts the observational range of cr = 30° -50°. The PP-i-FLRW 
model (red curve) can reproduce the observed SEP event widths, 
with width cr = 34°, while the FP model (blue curve) can not; 
it has a very narrow SEP width, with only cr - 11°. This strong 
difference between the FPh-FLRW and FP models is consistent 
with Laitinen et al. (2013|l, who found for Cartesian geometry 


that at early times the mean square width in full-orbit simula¬ 
tions exceeded the FP model result by an order of magnitude 
and that an FPh-FLRW model reprocudes the behaviour seen in 
full-orbit simulations well. 

Our results show the importance of modelling the early phase 
using a proper physical description. Both the FP and FPh-FLRW 
aim to describe the effect of the turbulence on the particles, in¬ 
cluding the cross-field spreading that is due to held-line me¬ 
andering. In the FP model, the field-line meandering is in¬ 
cluded only time-asymptotically using a diffusion description. 
However, as shown by Laitinen et al. (20131, the initial spread¬ 
ing of particles is fast and non-diffusive, relaxing to diffusive 


behaviour only at later times. This is accounted for in our new 
FPh-FLRW model and results in particles reaching much broader 
longitudinal extents than predicted by an FP model using the 
same turbulence parameters. 


Propagation along meandering held lines implies that similar 
event extents should also be reached with other particle species. 
This is indeed found in observations: |Lario et'af] ( |2013| ) and 
Richardson et al. (2014|l have found a similar range of peak in¬ 


tensity widths for electrons and protons at different energies, and 
[Wiedenbeck et al.| ( |2013| l and [Cohen et al.| ( [2014[ ) for heavy ions 
during impulsive events. While these observations do not ex¬ 
clude other SEP-spreading mechanisms, such as transport in the 
corona, symphatetic Haring and extended particle sources (see. 


e.g., Wiedenbeck et al. |2013[ Gomez-Herrero et al.j 2015 i, the 
effect of interplanetary Held-line meandering offers a direct ex¬ 
planation for the SEP event widths that are observed to be inde¬ 
pendent of energy and particle species and the gradual or impul¬ 
sive classiHcation. The influence of the source height, size, and 
different interplanetary conditions on the resulting SEP event ex¬ 
tent will be explored in a future full parametric study. We will 
also address the varying degrees of observed SEP anisotropies, 
which typically indicate cross-Held transport in the interplane¬ 
tary space (Dresing et al. 2014| ), but in some cases indicate rapid 
access of particles throughout the heliosphere ([Gomez-Herrero 
[etaLl[20T5l l. 
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Fig. 5. Contour plot of particle density in arbitrary units as a 
function of time and longitude at 1 AU radial distance from 
the Sun (red circle in Figs. 1 and 2). The top panel shows the 
FP-i-FLRW model result, while the bottom panel shows results 
for the FP model. 



Fig. 6. Peak intensity at 1 AU as a function of longitude for the 
FP (blue solid curve) and FP-hFLRW (red solid curve) models. 
The dashed lines show the fits of Gaussian profiles to the peak 
intensities, with cr = 11° for the FP model, and cr - 34° for the 
FP-i-FLRW model. The grey area depicts the range of Gaussian 
widths of 30 - 50° obtained from observations. 


4. Conclusions 

We developed a new SEP transport model that takes the non- 
diffusive propagation of SEPs early in the event history into ac¬ 
count for a Parker spiral geometry. We showed that the early on¬ 
set of SEPs over a wide range of longitudes can be explained by 
field-line random walk and requires an SEP transport model that 
properly describes the non-diffusive early phase of SEP cross- 
field propagation. Our EPh-EERW model is the first model that 


is capable of reproducing the observed fast access of SEPs to 
distant longitudes, when the particle and field-line diffusion co¬ 
efficients are consistently derived from an interplanetary turbu¬ 
lence model. When the ELRW is not included (in the EP model), 
a much narrower cross-field extent of the SEP event is pro¬ 
duced. We conclude that introducing field-line wandering into 
SEP modelling has the potential of resolving the problem of fast 
access of SEPs to a wide range of longitudes. 


Appendix A: Stochastic differential equations in the 
FP and FP+FLRW models 


The method of stochastic differential equations (SDEs) is a 
tool for solving Eokker-Planck-type transport equations. Instead 
of solving this equation directly, as is the case when apply¬ 
ing finite-difference methods, for example, the SDEs trace so- 
called pseudo-particles, that is, phase-space elements, through 
the phase-space. These pseudo-particles obey the corresponding 
Langevin equation. The simplest example of a transport equation 
is 


df df ^ 

dt dx dx^ 

(A.l) 

with the Langevin equation 


dx - V dt + yjlD dt W,, 

(A.2) 


where f(x, f) is the distribution function that depends on space 
X and time t. v and D denote velocity and diffusion coefficient, 
dt is the time step and and W, a Gaussian random number with 
zero mean and unit variance. The solution of the entire equa¬ 
tion is obtained by binning the results for an appropriate number 
of pseudo-particles. Although these phase-space elements repre¬ 
sent the behaviour of real particles in most cases quite well, they 
should not be confused with the latter. In contrast to the tradi¬ 
tional approach, the SDE method also allows solving the trans¬ 
port equation, with a different Langevin equation, backwards in 
time. Eor details on these points and on the method in general, 
we refer to Kopp et al. (20121 and the literature cited there. 

In the present case we solved the transport equation, Eq. Q 
under the assumption of a pure radial and constant solar wind 
(Vjw = Fj^Cr) and a vanishing latitudinal component of the mag¬ 
netic field, which in particular significantly reduces the Skilling 
term R (the term in curved brackets in the second line of Eq. Q 
). We first solved the equation in a local Cartesian system with 
z being the direction along the magnetic field, and transformed 
the transport quantities back to the global system before the ac¬ 
tual integration was carried out (Effenberger et al. 2012| l. In this 
local system, the Langevin equations read 


dx 


dy 

d^ 




dt 

dx 


dK_i_ y dt 

dy 


j t/f + ^t,x. 


1 

2 ’ 


i^-+^iRy\-^?)dt+ ^ dt + yiD^^dt W,,^.(A3) 


Here, u\\ - Vswb,-, where br is the radial component of the 
unit vector b = e^ along the local magnetic field line and 
u± = Vswi'^r - ^/b)- The two local unit vectors perpendicular to 
the magnetic field are x and e,. = x e^, and the corre¬ 

sponding particle cross-field diffusion coefficients » and a:j^ 
respectively. Eor the EP model, b is the unit vector along the lo¬ 
cal Parker field direction, whereas for the EP-i-PLRW model the 
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pseudo-particle is first propagated (virtually) along the mean¬ 
dering field line, and b itself is computed through the difference 
between the old and new points. 

The transport coefficients in Eq. ( |A.3[ ) are transformed into 
the global system by means of the tensor Tap - Cap, so that 
the (spatial) diffusion tensor becomes non-diagonal. In this case. 


the square root of a tensor has to be computed, see Kopp et al. 

dMTIl l. 

It is important to note that calculating the meandering field 
line using Eq. ([^ is performed before the particle is injected into 
the meandering field line. The particle cross-field propagation, 
while diffusive relative to the meandering field line, is therefore 
not diffusive relative to the mean magnetic field direction until 
the cross-field distance that is due to integration of Eq. ( |A.3| l is 
comparable to that due to the integration of Eq. (§. At early 
times, the perpendicular pseudo-particle propagation is there¬ 
fore only weakly stochastic with respect to the mean magnetic 
field, making our model consistent with the results presented by 
[Laitinen et al.| ( |^13[ |. 
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